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ABSTRACT 

We study about ultra-high-energy cosmic rays (UHECRs) from transient sources, propagating in the Galac- 
tic and intergalactic space. Based on the recent observational results, we also estimate upper and lower 
bounds on the rate of transient UHECR sources and required isotropic cosmic-ray energy input per burst as 
0. 1 Gpc"'' yr~' < po < lO'''^ Gpc"'' yr~' and lO^'''^ ergs < f hecr ^ ^i^gs, through constraining the apparent 
burst duration, i.e., dispersion in arrival times of UHECRs. Based on these bounds, we discuss implications for 
proposed candidates such as gamma-ray bursts and active galactic nuclei. 

Subject headings: cosmic rays — gamma rays: bursts — galaxies: active 



1. INTRODUCTION 

The origin of ultra-high-energy cosmic rays (UHECRs) is 
one of the biggest mysteries in astroparticle physics. So far, a 
number of possibilities were proposed, and several accelera- 
tion me chanisms have been theoretically developed (see, e.g., 
iKachelriefi 2008 , and references there in). However, physical 
conditions in these potential sources are uncertain, and obser- 
vational progress for source identification has been limited by 
the scarcit y of experimenta l data at the highest energies (see, 
e.g.. lNagano & Watson 2000i) . 

The recent observational results of large area detectors 
such as the Akeno Giant Air Shower Array (AGASA), 
High Resolution Fly's Eye (HiRes), and especially the 
Pierre Auger Southern Observatory (PAO), have started 
to give us crucial clues to the association of UHECRs 
with UHECR sources. Indeed, the first results of the 
PAO reported a significant correlation between the ar- 
rival directions of the highest-energy cosmic rays with 
the 12th Veron-Cetty & Veron catal og of active galactic 
nuclei (AGNs) closer than 75 Mpc (lAbraham et al. 20071: 
[Abraham et al. 2008) . Although thi s resul t has not been con- 
firmed by the HiRes (jAb basi et a l. 20081) and criticized by 
several authors dGorbunov et al. 2007h . it has also received 
some confirmations, and it would be an im portant step t o- 
wards solving the UHECR mystery (see, e.g. JStanev 20 08). 

However, one should not overinterpret the significance of 
these results. Although several authors also reported cor- 
relations of UHECRs w ith AGNs (iGeorge et al. 20081: 
iMoskalenko et al. 20081: IZaw et al. 2008h . one can- 
not exclude the possibility of other objects associated 
with the large scale structure of the universe, which 
is inhomogeneous up to dozens of Mpc. Significant 
correlatio ns of UHECRs with ga laxies can also be 
found dKashti & Waxman 20081: iGhiseUini et al. 20081: 
iTakami et al. 2008), so that gamma- ray burst s (GRBs) 
dVietri 1995; Waxma n 1995a; Mur ase et al. 2006 ) and 
magnetars (lArons 2003b can be sources. 

Even if the association of UHECRs with AGNs is real, the 
report by the PAO brought us several new questions on the 
nature of AGNs generating UHECRs. Surprisingly, the large 
majority of the correlating AGNs seems radio-quiet, a class 
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of objects not showing any nonthermal high-en ergy emission 
in their photon spectrum (iGeorge et al. 20081) . Radio-loud 
AGNs, showing high-energy nonthermal emission, are more 
plausible candidates in th e conventional jet paradigm (e.g., 
iRachen & Biermann 19931: iNorman et al . 1 9 950 . Althoug h 
the association with them is argued (Moskalen ko et al. 20081) . 
it seems that the power of th e correlatin g AGNs are insuf- 
ficient to produce UHECRs (IZaw et al. 2008) . The above 
problem may be solved if UHE CRs are produced durin g 
intense but short-duration flares (i Farrar & Gruzinov 20081) . 
The magnetic fields in the universe deflect UHECRs, 
so that UHECRs are significantly delayed compared 
to photons and neutrinos gene rated during the bursts 
jMiralda-Escude & Waxman 1996h . A transient hypothesis 
might also help to reproduce the isotro py of the arrival d is- 
ti-ibution of UHECRs at - 10''' eV (Tak ami & Sato 2008bh . 

In this letter, we focus on the possibility that UHECR 
sources are transient, and evaluate the deflection angles and 
arrival times of UHECRs through numerical calculations, 
considering both of the Galactic magnetic field (GMF) and 
intergalactic magnetic field (IGMF). The required cosmic-ray 
energy input and rate of the sources are estimated. In this 
work, UHECRs are also assumed to consist of protons. 

2. PROPAGATION AND CHARACTERISTICS OF UHECRS FROM 
TRANSIENT SOURCES 

We briefly describe the method of calculation and charac- 
teristics of UHECRs from transient sources. UHECRs ejected 
from their sources are deflected by the GMF and IGMF dur- 
ing their propagation. Only if the deflection angle dd{E,D) is 
small, where E is the energy of UHECRs at the Earth and D 
is the source distance, we could see a positional correlation of 
the highest-energy events with the sources at observationally 
suggested small-angle separations. The deflection also causes 
the time delay t^(E,D) between arriving times of an UHECR 
and a light emitted at the same time. UHECRs with the same 
energy have different arrival times, because of not only differ- 
ent parti cle trajectories but also stochast ic photomeson pro- 
duction (iMiralda-Escude & Waxman 199 6). Therefore, the 
time delay has certain distribution with an averaged delayed 
time td{E.D) and standard deviation in arrival times ad(E,D), 
The arrival time spread a-j can be regarded as the apparent 
duration of an UHECR burst. 

Clearly, the magnetic fields play an essential role on both 
of 9d and crj. For intergalactic propagation, one can typi- 
cally expect (Trf ~ frf « ^ ~ 10^ yrs £'20^iooMpc'BiG,-9'^Mpc 
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dMiralda-Escude & Waxman 19961) . which is also confirmed 
by our numerical calculations. Due to limited statistics of 
the highest-energy events, it is convenient to use quantities 
weighted by the observed cosmic -ray spectrum. The apparent 
burst duration of UHECRs above the threshold energy Eth is 



Td(>Eih)-- 



1 



3 J A r r^" 



x(£) 
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dE 
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(1) 



where dNo/dE is the UHECR spectrum observed at the 
Earth, A/q = dE^^{E) is the normalization factor, and 
DmaxiE) is the maximum distance of UHECRs that can reach 
the Earth at the energy E. In this work, we adopt E± = lO'^'^^ 
eV as the threshold energy, according to the R^O results. 

Through t^, we can relate the local rate of transient sources 
po with the apparent source density n.,. We have 



n,(> -Eth) ■■ 



1 



Eih 
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where 
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(3) 



Note that can be estimated from the observed small 
scale anisotropy in the arrival distribution of the highest- 
energy cosmic rays with energies above E^^. For example, 
the small scale anisotropy observed by the AGASA im- 
phed 10-''-10-^ Mpc" ^ (e.g., lYoshiguch i et al. 2003t 

iKachelrieB & Semikoz 200l |j akami & Sato 20071) . 
The more recent PA O data imply n, ^ 10"^ Mpc"'' 
dTakami & S ato 2008bl). and we hereafter adopt this value. 
Then, the local burst rate is estimated via po ~ ris/T^. 

Assuming that the sources are uniform, we can also esti- 
mate typical values of the isotropic cosmic-ray energy input 

per burst at the energy E as £'(^^{E) w E^^{E)/pa. Here, 

E^^2^(E) is the UHECR energy budget per volume per year 
at the energy E. Through our numerical calculations, we ob- 
tain £2!i^(ioi9 eV) ~ (0.5-2) x 10*"* ergs Mpc"^ yr"' (de- 
pending on the source spectral index s) from the PAO data 
(see also, e.g., Waxman 1995b; Berezinskv et al. 2006). 

The thing left to do is to calculate the distribution of deflec- 
tion angles of arrival times. Our method of calculation, taking 
into account the GMF as well as the IGMF, is described as 
below. The IGMF strength Big is very uncertain, but we can 
estimate upper bounds on Big and resulting r^, by comparing 
the calculated distribution of deflection angles to the typical 
angular separation of observed UHECRs. In this work, we 
adopt ^ ~ 5° as the angular separation, according to the PAO 
results (Abraha m et al. 20081: fTakami et al. 20081) . As the en- 
ergy distribution of cosmic rays at a source, we assume power- 
law spectra of dN/dEg oc E'". 

2.1. Propagation in the Galactic Space 

Propagation in the Galactic space is important for the de- 
flection of UHECRs. The corresponding delay time would 
typically be smaller than that by the IGMF, but it is not zero 
and is unavoidable. Under a given separation angle, the co- 
herent component of the GMF leads to the minimum delay 
time, and the lower bound on is also obtained. Following 
the method used in Takami & Sato 2008a, we pursue cosmic- 
ray trajectories with proton mass and charge of -1 from the 



TABLE 1 

Upper bounds on t^' from the GMF without a dipole field. 



.V 


-r,7' [yr'] forp = -10° 


r;' [yr'] forp = -8° 


2.0 


2.5 X 10-- 


3.1 X 10-2 


2.2 


2.5 X 10"- 


3.0 X 10-2 


2.4 


2.5 X 10-- 


3.0 X 10-2 


2.6 


2.4 X 10-- 


2.9 X 10-2 



Upper bounds on ' 



TABLE 2 

FROM THE GMF WITH A DIPOLE FIELD. 



.V 


Tj' [yi-'] forp = -10° 


rji [yr-']forp = -8° 


2.0 


6.4 X 10^ 


6.1 X 10-"* 


2.2 


6.3 X K)-^ 


6.0 X 10-* 


2.4 


6.2 X 10^ 


5.9 X 10^ 


2.6 


6.2 X 10^ 


5.9 X 10~* 



Earth, and calculate their delay times for given GMF mod- 
els. We define a sphere with the radius of 40 kpc, centered of 
the Galactic center, as the boundary of Galac tic space. As a 
GMF model (see, reviews, IVallee 20041: iHan 2007i) . a bisym- 
metric spiral field with the even parity is adopted. This leads 
to conservative estimate of t^, although other models such as 
axisymmetric spiral field may be possible due to uncertainty 
in the GMF ( Vallee 2005). Note that all the energy loss pro- 
cesses can be neglected for propagation in the Galactic space. 

The delay time and arrival time spread by the GMF depend 
on arrival directions of UHECRs, where the averaged stan- 
dard deviation of observed k events is crj = (1/A:)S,(t,/., . In this 
work, we instead use (7d(E) = J dfl ^(E,n). Through our 
numerical calculations, we found that this gives us reasonable 
estimate of (t^, even though the GMF leads to the h ole in ar- 
rival directions of UHECRs (iTakami & Sato 2008ah . In Table 
1, we show resulting lower bounds on tj. Here p (< 0) is the 
pitch angle of the spiral component of the GMF at the vicinity 
of the solar system, and smaller values of -p leads to smaller 
deflection angles of UHECRs. 

The vertical magnetic field near the solar system and many 
gaseous filaments perpendicular to the Galactic plane are 
observed, indicating another regular component (Han 2007). 
Also, a dipole field with the odd party is predicted by the dy- 
namo theory. Hence, we also consider cases of the GMF with 
a dipole magnetic field whose strength is normalized to 0.3 
fiG at the vicinity of the solar system. In Table 2, resulting 
lower bounds on tj are shown. However, note that there is no 
direct observational evidence of such a dipole field. 

2.2. Propagation in the Extragalactic Space 

Propagation in the extragalactic space is expected to play 
an essential role on both of the deflection and time delay of 
UHECRs. We numerically calculate the distribution of and 
td for given IGMFs. Our method of calculation is similar to 
that used in Yoshiguchi et al. 2003, where proton propagation 
is treated as Monte Carlo simulations. We set 10 logarithmic 
bins per logarithmic energy interval, and isotropically inject 
5000 protons for every energy bin. Particle trajectories are 
pursued at every 1 Mpc for D < 100 Mpc while at every 10 
Mpc for D > 100 Mpc. As relevant energy loss processes, 
we consider photomeson production and Bethe-Heitler pro- 
cesses with the cosmic microwave background photons, and 
the adiabatic energy loss due to the expanding universe. The 
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TABLE 3 

Obtained values of t^' for uniform turbulence IGMFs. 



s 


-r,;' [yr'] for 0.1 nGMpcl/2 


r;' [yr'] for l.OnGMpc'/^ 


2.0 


9.7 X 10-5 


8.9 X 10-' 


2.2 


9.5 X 10-^ 


8.6 X 10-' 


2.4 


9.3 X 10-'' 


8.4 X 10-' 


2.6 


9.1 X 10-^ 


8.3 X 10-' 



resulting distribution of 64 or is 



dEg fd(E,D;Eg), 



(4) 



where £'™"(£'th,£') is the minimum energy at a source with 
the distance D, of protons observed with Ei^ at the Earth. 
fd{E ,D;Eg) is the more basic distribution of 6^ or f^/, gen- 
erated by cosmic rays with Eg at a source. 

In this work, we consider an uniform turbulent IGMF with 
the Kolomogorov turbulence spectrum as the extragalactic 
magnetic field. Although the IGMF is highly uncertain, we 
can constrain it by comparing the calculated 0^ distribution 

to ip. As a result, we found BigA'^^ < nGMpc'^^ is required 
for the averaged deflection angle 9d not to exceed the typi- 
cal angular separation ijj 5°. It is consis tent with the re - 
sult from Faraday rotation measurements (Kronberg 19941) . 

1/2 

but independently obtained. Hence, obtained for Big Aj,^jj ~ 
nGMpc'/^ can be regarded as upper bounds. Our results for 
Big = (0.1 - 1) nG and Acoh = 1 Mpc are shown in Table 3. 

3. IMPLICATIONS FOR TRANSIENT UHECR SOURCES 

We have estimated lower and upper bounds on using ip ^ 
5°, which allows us to estimate the allowed range of po and 
^CR using ris ^ 10""^ Mpc""', As the local rate, we obtain 



0.1Gpc"'yr- 



"' <Po<(60- 



3000)Gpc-'yr- 



Note that stronger upper bounds can be obtained with the 
GMF with a dipole magnetic field. However, since the exis- 
tence of a dipole field is very tentative, we hereafter consider 
the GMF without a dipole field for conservative discussions. 



The required cosmic -ray energy input at lO'^ eV, ^hecr 



F' 



Pjj(10'^ eV) is estimated as 
(0.3 



20) X 10^° ergs < 4ecr ^ 10^"* ergs. 



(6) 



Note that Eqs. (5) and (6) are valid as long as UHECR sources 
are regarded as transient, i.e., 6T <Td < AT, where 6T is the 
true burst duration during which particle acceleration occurs 
and Ar is the time interval between bursts. ST depends on 
the nature of potential sources. For example, classical high- 
luminosity (HL) GRBs have ST ^ 10'~^ s, which is much 
shorter than rj. For bounds to be meaningful, Td < AT should 
be satisfied. Otherwise, more than one UHECR bursts oc- 
cur in Td within ip, and we would see these bursts as a single 
but more energetic burst. When ^ ^ 5°, the time interval is 
estimated as Ar ~ (3/7r)pQ V~^^^max(^th) — 3t^«71_4- Since 
Ar > Td, we may expect that obtained bounds would make a 
sense, although we should be careful of the possibility not to 
see each UHECR burst as a distinctive one for larger n,. 

The total cosmic-ray energy input is generally larger 
thang-cRbyj^dO '^ eV) ^ (J dE'g E'g / (E^ §-) ^^^^^^^ 
dMurase et al. 20081) . R depends on the cosmic-ray spectrum 



TABLE 4 

Potential Sources and Their Typical Local Rates 



Source 


' — 1 1 — 

Typical Rate po [Gpc yr ] 


Reference 


HLGRB 


~0.1 


e.g., GP07 


LLGRB 


~400 


e.g., L+07 


Hypemovae 


~ 2000 


e.g., GD07 


Magnetar 


~ 12000 


e.g., G+05 


Giant Magnetar Flare 


~ 10000 


e.g., O07 


Giant AGN Flare 


~ 1000 


FG08 


SNe Ibc 


~ 20000 


e.g., GD07 


Core Collapse SNe 


120000 


e.g., M+98 



at a source, and we expect R ~ 20-500 for s 2.0-2.2 ex- 
pected in the ankle scenario while 7? > 100 for s ^ 2.4-2.6 
expected in the proton-dip scenario, and the latter scenario 
generally requires the break energy below the second knee 
jBerezinskv et al. 200 6). In both scenarios, we expect that 
the transient hypothesis requires the relatively large cosmic- 
ray energy input per burst > 10^"'^ ergs, which would be 
a strong requirement on potential sources. 

So far, several potential sources are proposed as transient 
accelerators, and HL GRB is one of them. The isotropic 



radiation energy is ^ 10^^ ergs, and HL GRBs are the 
most energetic transient phenomena in the universe. The 
local rate is uncertain, but recently suggested rates in the 
Swift era, po (0.05-0.27) Gp c ~^ yr~' are smaller t han 
previous ones jLe & Dermer 20071: iGuetta & Piran 2007h . If 
Po < 0.1 Gpc"^ yr-> is real, HL GRBs would be diffi- 
cult as UHECR sources, since they require rather strong 
IGMFs with Big ^ nG and large isotropic energy in- 
put of f > 2 X 1Q^\R/20) ergs. Low-luminosity (LL) 
gamma-ray bursts may overcome the p roblem that the lo- 
cal rate of HL GR Bs seems too small (iMurase et al. 20061 : 
iMurase et al. 20081) . because their local rate is Ukely to 
(5) be much higher, pp ^ 10^ "^ Gpc"^ yr"' ('Liang et al. 20071: 



iGuetta & Delia Valle 200% . Some GRBs ai-e associated with 
energetic supernovae (SNe) called hypernovae, which may 
also be high-energy cosmic-ray accelerators. Although they 
are sufficient as the energy budget, it seems difficult to accel- 
erate protons up to > 10'^ eV dWang et al. 200 7). 

Ab out 10 % of core co llapse SNe may form magnetars 
(e.g., iGaensler et al. 2005h . which may be UHECR sources 
dArons 2003 f). However, our results would suggest that all 
the magnetars (and SNe) do not produce UHECRs uniformly 
and only a fraction of magnetars is the main origi n, which i s 
also consistent with the theoretical expectation ( A rons 2003h . 
For example, only newly born magnetars associated with 
SNe Ibc could be major UHECR accelerators, which leads 
to po - 3000 Gpc"^ yr-i and f - 3 x 10^"(7;/10) ergs. Gi- 
ant magnetar flares could not explain UHECRs, since their 
radiation energy, £1!^° 10"*^ ergs is much smaller than 

AGNs are the most discussed UHECR accelerators, and 
the possibility as transient sources was recently suggested 
(iFarrar & Gruzinov 2008). Although such giant AGN flares 
may be UHECR sources, the suggestion is speculative. 
If we adopt po ~ 10^"^ Gpc~^ yr"' although the rate is 
also uncertain, the required energy input is S^^ ~ 2 x 



2q51 52|-^^20) ergs. The corresponding luminosity is L^^ 



2 X 10^""^' {R/20) (10' s/6T) ergs s"'. 

The typical source rates are summarized in Table 4, but one 
should keep in mind that they are very uncertain at present. 
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4. SUMMARY AND DISCUSSIONS 

We have constrained through our numetical calculations 
of proton propagation, taking into account not only the IGMF 
but also GMF. The recent PAO results, the positional correla- 
tion with a few degrees and small-scale anisotropy, have 
brought us implications for the sources, as well as indicated 

Z?igA^^^ < nG Mpc'''^. The results would be important in 
the sense that they are implications of UHECR observations 
and also useful for the detectability of secondary emission. 
They suggest that HL GRBs may be marginally disfavored as 
UHECR sources, if the recently suggested local rate is real. 
LL GRBs, newly born magnetars associated with LL GRBs 
or SNe Ibc, and giant AGN flares seems possible, unless there 
exists a strong dipole magnetic field in our Galaxy. 

Although these implications may be interesting, because of 
current poor stastistics, we should be careful to reach defi- 
nite conclusions about the sources. The suggested posional 
correlation is confirmed just at 2-3 sigma levels. Also, the es- 
timate of «s has large errors at present due to not only poor 
statistics but also the lack of our knowledge of the precise po- 
sitions of UHECR sources, and it is not so easy to exclude the 
possibility that UHECR sources contains many dim acceler- 
ators, i.e., large n^. But, statistics will be better in the near 
future, and the order of will be accurately determ ined by 
5 yrs observations by the PAO jTakami & Sato 2007h . How- 
ever, statistical analyses become more complicated when the 
sources are transient and/or their luminosity function is taken 
into account (and note that we have assumed that the sources 
are uniform). More detailed and careful studies will be pre- 
sented in our forthcoming paper. 

We have discussed implications for the transient UHECR 
sources, but it is also important to know whether the sources 
are transient or not. The signature of transient sources may 
be found from the observed UHE CR spectrum, e.g. , from 
the average number of multiple ts dHarari et al. 2004| ). For 
identification of sources, the multimessenger astronomy will 
be particularly important. High-energy neutrinos are use- 
ful as a probe of cosmic -ray acceleration, and associated 



gamma rays are oft en expected f o r transient sources such 
as GRBs (see, e.g., iMurase 2007t iMurase et al. 20081 and 
references there in). In addition, electrons are also accel- 
erated in the shock acceleration theory, which allows us to 
expect photon counterparts. In fact, we may also expect 
the high electi-on luminosity L, > (e,/0.1ep)10'^^(10^ s/6T) 
ergs s"', as well as the high magnetic luminosity Lb > 
lO"^^ ergs s~^ whi ch will be required for UHECR acceleration 
jWaxman I995al) . For example, ~ 30 events from giant AGN 
flares may be detected by Fermi (iFarrar & Gruzinov20 08). 
Note that, once we know that the UHECR sources are tran- 
sient, w e can obtain precious informati on on the effective 
IGMF (Miralda- Escude & Waxman 1996 ). If LL GRBs are 
the UHECR sources, for example, the effective IGMF can be 
estimated as - 0.03 nG Mpc'/^ Also, UHECRs would also 
be useful as a probe of the GMF (T akami & Sato 2008al) . and 
several GMF models may be tested as the number of detected 
events increases. This can help to reduce uncertainties in es- 
timate of bounds on and resulting po. 

We have assumed an uniform IGMF, but it would not be 
realistic. If the sources are inside clusters or filaments, UHE- 
CRs are affected by the magnetic field in the structured region. 
When an uniform IGMF is weak enough, such the structured 
magnetic field is more important. In such cases, HL GRBs are 
more disfavored as the UHECR sources because they require 
the rather strong IGMF, as shown in this work. The structured 
magnetic field may also play a role as an unavoidable field. 
Some AGNs are inside clusters, so that UHECRs from them 
should be delayed due to that field. The local rate will be 
more constrained as po < 10'~^ Gpc~^ yr"', and the necessary 
cosmic-ray energy input will also be increased. 
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